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^ ' Abstract 

oo : 

ly-) t We obtained steady solutions of optically thin, single temperature, magnetized black hole accretion disks 

\Q . assuming thermal bremsstrahlung cooling. Based on the results of 3D MHD simulations of accretion disks, 

we assumed that the magnetic fields inside the disk are turbulent and dominated by azimuthal component. 
We decomposed magnetic fields into an azimuthally averaged mean field and fluctuating fields. We also 
assumed that the azimuthally averaged Maxwell stress is proportional to the total pressure. The radial 
advection rate of the azimuthal magnetic flux $ is prescribed as being proportional to vj~^ , where zu is 
■ the radial coordinate and £ is a parameter which parameterizes the radial variation of 4>. We found that 

when accretion rate M exceeds the threshold for the onset of the thermal instability, a magnetic pressure 
q ■ dominated new branch appears. Thus the thermal equilibrium curve of optically thin disk has a 'Z'-shape 

\ in the plane of surface density and temperature. This indicates that as the mass accretion rate increases, 

c/3 . a gas pressure dominated optically thin hot accretion disk undergoes a transition to a magnetic pressure 

dominated, optically thin cool disk. This disk corresponds to the X-ray hard, luminous disk in black hole 
, candidates observed during the transition from a low/hard state to a high/soft state. We also obtained 

global steady transonic solutions containing such a transition layer. 
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1. Introduction flows (RIAFs), were studied extensively by Narayan, Yi 

(1994,1995) and Abramowicz et al. (1995). Abramowicz et 

Optically thin, hot accretion disks have been stud- al. (1995) obtained local thermal equilibrium curves of op- 

ied to explain the X-ray hard state (low/hard state) of tically thin accretion disks by solving Q v i s = Q a dv + <3rad, 

black hole candidates. Thorne, Price (1975) proposed where Q rac j is the radiative cooling rate. They showed 

that hard X-rays from Cyg X-l are produced in the in- that when the accretion rate exceeds a critical accretion 

ncr optically thin hot disks. Shibazaki, Hoshi (1975) rate, optically thin solution disappears. Above this criti- 

studied the structure and stability of optically thin hot cal accretion rate, the hard X-ray emitting optically thin 

accretion disks. Eardley, Lightman, Shapiro (1975) and hot accretion disk undergoes a transition to an optically 

Shapiro, Lightman, Eardley (1976) constructed a model thick cool disk which emits soft X-rays. The latter corre- 

of optically thin, two temperature accretion disks in which sponds to the high/soft state (or thermal state) of black 

ion temperature is higher than the electron tempera- hole candidates. 

ture. Ichimaru (1977) pointed out the importance of The importance of the magnetic fields on transition be- 

advection of energy in hot accretion flows and obtained tween an optically thin hot disk and an optically thick 

steady models of optically thin disks by equating the ra- cool disk was discussed by Ichimaru (1977). He sug- 

dial heat advection Q a dv and the viscous heating Q v i s . gested that magnetic pressure p m ag is limited below the 

Such flows, which now arc called advection dominated ac- gas pressure p gas because magnetic flux will escape from 

cretion flows (ADAFs) or radiatively inefficient accretion the disk by buoyancy. However, the buoyant escape of the 
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magnetic flux can be suppressed in strongly magnetized 
disks. Shibata, Tajima, Matsumoto (1990) carried out 
magnctohydrodynamic (MHD) simulations of the buoy- 
ant escape of the magnetic flux due to the Parker in- 
stability (Parker 1966) and showed that when the disk 
is dominated by magnetic pressure, it can stay in low-/? 
state (/3 = Pgas/Pmag < 1) because the growth rate of the 
Parker instability is reduced in low-/3 disks due to the mag- 
netic tension. Mineshige, Kusunose, Matsumoto (1995) 
suggested that low-/3 disk emits hard X-rays. Pariev, 
Blackman, Boldyrcv (2003) obtained steady solutions of 
optically thick, magnetic pressure dominated disk but 
they did not consider the optically thin solution including 
advective cooling. Begelman, Pringle (2006) constructed 
a model of optically thick low-/? disk including radiation 
pressure. 

Since magnetic fields contribute both to the angular 
momentum transport and disk heating, we need to under- 
stand how magnetic fields are amplified and maintained 
in the disk. In conventional theory of accretion disks (e.g., 
Shakura, Sunyaev 1973), phcnomenological a- viscosity is 
invoked. However, the physical mechanism which en- 
ables angular momentum transport efficient enough to 
explain the activities of accretion powered sources such 
as dwarf nova was unknown until Balbus, Hawley (1991) 
pointed out the importance of the magneto-rotational in- 
stability (MRI) in accretion disks. MRI can excite and 
maintain magnetic turbulence. The Maxwell stress gen- 
crated by the turbulent magnetic fields efficiently trans- 
ports angular momentum and enables accretion of the 
disk material. The growth and the saturation of MRI 
has been studied by local three-dimensional MHD simula- 
tions (e.g., Hawley et al. 1995, 1996; Matsumoto, Tajima 
1995; Brandenburg et al. 1995; Sano, Inutuska 2001; Sano 
et al. 2004) and by global three-dimensional MHD simula- 
tions (e.g., Matsumoto 1999; Hawley 2000; Machida et al. 
2000; Hawley, Krolik 2001; Hawley, Krolik 2002; Machida, 
Matsumoto 2003). These global three-dimensional MHD 
simulations of radiativcly inefficient accretion disks indi- 
cate that the amplification of magnetic fields saturates 
when /? rsj 10 except in the plunging region of black hole 
accretion disk and in the disk corona, and that the disk 
approaches to a quasi-steady state. In this quasi-steady 
state, the average ratio of the Maxwell stress to gas 
pressure ass, which corresponds to the a-parameter in 
the conventional accretion disk theory (Shakura, Sunyaev 
1973), is - 0.01-0.1 (e.g., Hawley 2000; Hawley, Krolik 
2001; Machida, Matsumoto 2003). 

RXTE observations of state transitions in galactic black 
hole candidates revealed that some black hole candidates 
(e.g., GX339-4, XTE J1550-564, and XTE J1859+226) 
stay in X-ray hard states even when their luminosities 
exceed 20% of the Eddington luminosity (e.g., Done, 
Gierlihski 2003; Gierlihski, Newton 2006). When a < 0.1 
as indicated by global three-dimensional MHD simulations 
of black hole accretion flows and the energy conversion ef- 
ficiency r) c ~ 0.1, their luminosities are higher than the 
critical luminosity above which optically thin steady hot 
solutions disappear. 



Recently, Machida, Nakamura, Matsumoto (2006) car- 
ried out global three-dimensional MHD simulations of 
black hole accretion disks including optically thin radia- 
tive cooling. They started the simulation from a radia- 
tivcly inefficient torus threaded by weak toroidal mag- 
netic fields. As MRI grows, an optically thin, hot accre- 
tion disk is formed by transporting angular momentum 
efficiently They found by simulation including thermal 
brcmsstrahlung cooling that when the density of the ac- 
cretion disk exceeds the critical density, cooling instabil- 
ity takes place in the disk. The disk shrinks in the verti- 
cal direction with almost conserving the toroidal magnetic 
flux. They demonstrated that as the disk shrinks, mag- 
netic pressure exceeds the gas pressure because the mag- 
netic pressure increases due to flux conservation, mean- 
while the gas pressure decreases due to cooling. When 
the magnetic pressure becomes dominant, the disk stops 
shrinking in the vertical direction because the magnetic 
pressure supports the disk. They found that the disk 
evolves toward a quasi-equilibrium cool state. During this 
transition from hot state to cool state, the disk remains 
optically thin. They explained this quasi-steady state by 
equating the magnetic heating rate Q + and the radia- 
tive cooing rate Q la d assuming that the Maxwell stress is 
proportional to the total pressure and that the azimuthal 
magnetic flux (B V )H is conserved during the transition, 
where (B v ) is the mean azimuthal magnetic field, and H 
is the half thickness of the disk. 

The purpose of this paper is to answer why hard X-ray 
emitting optically thin disks exist above the critical lumi- 
nosity. We show steady solutions of optically thin black 
hole accretion disks by taking into account the magnetic 
fields, radiative cooling, and advection of energy and mag- 
netic fields. In section 2, we present basic equations. In 
section 3, we obtain local thermal equilibrium curves. We 
present the global steady solutions in section 4. Section 5 
is devoted for summary and discussion. 

2. Models and Assumptions 

2.1. Basic Equations 

We extended the basic equations for one-dimensional 
steady, optically thin black hole accretion flows (e.g., Kato 
et al. 1998) by incorporating the magnetic fields. We 
adopt cylindrical coordinates (zu,ip,z). General relativis- 
tic effects arc simulated by using the pseudo-Newtonian 
potential tp = —GM/(r — r s ) (Paczynsky, Wiita 1980), 
where G is the gravitational constant, M is the mass 
of the black hole (we assume M = 10M Q in this paper), 
r = (w 2 + z 2 ) 1/2 , and r s = 2GM/c 2 is the Schwarzschild 
radius. 

We start from the resistive MHD equations 
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where p is the density, v is the velocity, B is the magnetic 
field, j = cV x J5/47T is the current density, p gas = 3?pT/ p 
is the gas pressure, T is the temperature, e = 3!ftT/2/i is 
the internal gas energy, 5ft is the gas constant, p is the 
mean molecular weight (we assumed to be 0.617), q + is 
the heating rate, q~ is the radiative cooling rate, and r\ is 
the resistivity. 

Global three-dimensional MHD simulations (e.g., 
Matsumoto 1999, Hawley 2000, Machida et al. 2000; 
Hawley, Krolik 2001, 2002; Machida, Matsumoto 2003; 
Machida et al. 2004, 2006 ) showed that in radiatively in- 
efficient accretion disks the amplified turbulent magnetic 
fields saturates when (3 = p gas /Pmag ~ 10 except in the 
plunging region and in the disk corona. Inside the disk, 
the azimuthal component of magnetic field dominates the 
poloidal component. 

Machida, Nakamura, Matsumoto (2006) carried out 
three-dimensional MHD simulation of black hole accretion 
disks including optically thin radiative cooling and showed 
that when the mass accretion rate exceeds the threshold 
for the onset of the cooling instability, a magnetically sup- 
ported disk is formed. The magnetic pressure exceeds the 
gas pressure because magnetic fields are amplified due to 
the vertical contraction of the disk, meanwhile the gas 
pressure decreases due to cooling. 

We have to take into account magnetic fields to study 
the evolution of accretion disks in this regime. We decom- 
pose the magnetic fields into the axisymmetric toroidal 
component (mean field) B = (B^e^ and the fluctuat- 



ing fields SB = 5B^e v 



SB v e v 



5B z e z and decompose 



the velocity into the mean velocity v = (u CT ,w v ,w 2 ) and 
fluctuating velocity Sv. Here ( ) denotes the azimuthal 



average and £? ro = SB^, B v = (B v ) 



5B Lpi and B z = SB Z 



We assume that the fluctuating components vanish when 
azimuthally averaged {(Sv) = (5B) = 0). Thus, the az- 
imuthally averaged magnetic field and velocity arc (B) = 
(B v )e v and (v) = v, respectively. 

Let us derive azimuthally averaged equations assuming 
that the disk is in a steady state. We assume that | (B v ) + 
SB V \ ^> \SB m \, \SB Z \. By azimuthally averaging equations 
(1) - (3) and ignoring the second order term of SB^ and 
SB Z , we obtain 
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[(pt+P gaB )v z ] 
= q + -q- ,(9) 



where p tot = p g 



dw dz 
- ftnag is the total pressure and p ma g 



(B^) /8ir is the azimuthally averaged magnetic pressure. 

Global three-dimensional MHD simulations showed 
that, in radiatively inefficient disks, azimuthally aver- 
aged total stress inside the disk is dominated by the wip- 
component of the Maxwell stress ~ (B^B V / Att) and 

that the average ratio of the Maxwell stress to the mag- 
netic pressure, a m = —(B^B v /4ir)/(p mlig ), is nearly con- 
stant (-0.2-0.5 ; e.g., Hawley, Krolik 2001, 2002; Pessah 
et al. 2006). Since /3 ~ 10 inside the disk except in the 
plunging region, a B = -(B^B^/4-k) /p to t ~ 0.01-0.1 (e.g., 
Hawley, Krolik 2001). Machida, Nakamura, Matsumoto 
(2006) showed that «b ~ 0.05 — 0.1 in cooling dominated 
magnetic pressure supported disk. According to their nu- 
merical results, we assume that the azimuthally averaged 
wip component of the Maxwell stress is proportional to the 
total pressure ((B^B v /4ir) = — asPtot)- Thus, equation 
(7) becomes 
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We note that in the turbulent accretion disk, the av- 
erage of the radial component of magnetic fields, (B m ), 
is small because positive and negative cancel out. 
On the other hand, the product of the radial and the 
toroidal components, does not change sign in the 

non- linear state (see figure 1). Thus, we can not neglect 
(B^B V ) term. 

2.2. Heating and Cooling Rates 



q + 



The other key factor is the heating rate expressed as 
t-^ipW (dfl j dw) in conventional theory, where fi is 
the angular velocity. Meanwhile, three-dimensional MHD 
simulations indicate that the total dissipative heating rate 
(a major part of dissipation is due to the thermalization 
of magnetic energy via magnetic rcconncction) is q + — 
(B^B ip /A'!T)w(dfl/dw) and that this dissipation occurs 
throughout the disk (e.g., Hirose et al. 2006). 

We assume the magnetic heating as the heating mech- 
anism in the disk and set the heating term as follows: 



(Brr, B,, 



dn 



dn 



-w— = -a B p tot w— . (11) 
47T dw dw 

As the cooling mechanism, we assume the radiative 
cooling by the optically thin thermal bremsstrahlung 
emission. Hence, the cooling rate is expressed as, 

q~ = 6.2 x 10 20 p 2 T 1/2 ergs s" 1 cm" 3 . (12) 
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2.3. Vertically Integrated Equations 

We assume that the temperature T, [3 (= p gas /Pmag), 
radial velocity and specific angular momentum £(= 
wv v ) are independent of z. Machida et al. (2000) and 
Miller, Stone (2000) showed that the magnetic loops 
emerge from the disk and form coronal magnetic loops. 
In this paper, we focus on the disk and ignored the low 
density and low (3 corona. We assume hydrostatic balance 
in the vertical direction and therefore ignore the left-hand 
side of the z-component of momentum equation. Under 
these assumptions, the surface density E, the vertically 
integrated pressure W, and the half thickness of the disk 
H are given by 



(13) 



pdz = v 2-k pq(vj)H 



r°° sqj 1 

w= / Ptot dz = — (l + zr 1 ^ 

J — oo A* 



H 



1 

where po(w) 



1/2 



1 



1/2 



(14) 



(15) 



is the equatorial density and f^K = 
(GMjw) 1 ! 2 j(vo — r s ) is the Keplerian angular speed. 

We now integrate other basic equations in the verti- 
cal direction. Since the density vanishes at the disk sur- 
face, the vertically integrated equation of continuity is ex- 
pressed as 



M 



(16) 



where M is the accretion rate. 

The ro-componcnt of the vertically integrated mo- 
mentum equation can be obtained by using (B v 2 ) = 

= V*/ (i + /r !) , 



dv 



1 dW 

dvj £ dzu 



Wdln^K 2 p- 1 W 



(17) 



m 3 £ dm w 1 + ft 1 £ 

where = n7 2 f2i<; is the Keplerian angular momentum 
and the second term on the right-hand side is a correction 
resulting from the fact that the radial component of the 
gravitational force changes with height ( see Matsumoto 
et al. 1984). 

Assuming that the last term of the (^-component of 
momentum equation vanishes at the disk surface, the (p- 
component of vertically integrated momentum equation 
is 

M{i-i- m ) = 2ira B m 2 W , (18) 

where l- m is the specific angular momentum swallowed by 
the black hole. 

The vertically integrated energy equation is 

Qadv = Q + -Q~ , (19) 
where 

M 9?T 



Q 



adv 



27TZZ7 2 fl 



< 



(20) 



is the advective cooling rate, where 



_ 3 d\nT + (ME dlnH 
2 dlnvo dlntu rflntu 



and 



q^dz = — aRWm— — , 
dw 



/oo 
q-dz = 6.2 x W 20 plT 1/2 y/^H 
-OO 



(21) 



(22) 



(23) 



are the heating and cooling rates integrated in the vertical 
direction. 

2-4- Prescription of the Advection Rate of the Toroidal 
Magnetic Flux 

The basic equations of the steady, magnetized accretion 
disks can be closed by specifying the radial distribution 
of magnetic field. By azimuthally averaging equation (4), 
we obtain 



d(B v )e v 
dt ~~ 



V x (v x (B v )e v ) 

+ Vx(5vxSB)+Vx (rjV x {B v )e v ) , (24) 

where the second term on the right-hand side is the dy- 
namo term and the last term is the magnetic diffusion 
term. If we neglect these terms and assume that the disk 
is in a steady state, we have 



(25) 



Assuming that the azimuthally averaged toroidal mag- 
netic fields vanish at the disk surface, we can integrate 
equation (25) in the vertical direction and obtain, 



v^(B v )dz = constant 



(26) 



The radial advection rate of the toroidal magnetic flux 
$ (hereafter we call it flux advection rate) can be evalu- 
ated by completing the vertical integration as 



(27) 



$ = -v^B (w)w AttH 
where 



B (m) = (B tp )(m;z = 0) 

= 2 5 /V/ 4 (9?T/ M ) 1 / 2 £ 1 / 2 ff- 1 / 2 /r 1/2 (28) 

is the azimuthally averaged equatorial toroidal magnetic 
field. In accretion disks, $ is not always conserved be- 
cause $ can change with radius by the presence of the 
dynamo term and the magnetic diffusion term in equa- 
tion (24). According to Machida, Nakamura, Matsumoto 
(2006), $ oc w^ 1 in the quasi steady state. Based on this 
simulation result, we parametrize the dependence of $ on 
w by introducing a parameter £ as 



( —) 



(29) 
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where ^out is the flux advection rate at the outer bound- 
ary w = tUout- When £ = 0, the magnetic flux is conserved. 
When £ > 0, the magnetic flux increases with approaching 
to the black hole. 

3. Local Model 

3.1. Local Approximation of Energy Equation 

Before obtaining global transonic solutions, we solve the 
energy equation Q + = Q ac iv + Q~ locally at some specified 
radius to obtain thermal equilibrium curves of optically 
thin black hole accretion flows. 

We approximated the energy equation in order to be 
solved locally. We evaluated £ by using results of global 
three-dimensional MHD simulations of black hole accre- 
tion flows. According to the results by Machida et al. 
(2004), T tx w' 1 , £ oc vo 1 ' 2 , and H cx w, thus £ = 1. 
Furthermore, we assume fl = fix- In our local model, we 
adopt £i n = ^K(3r s ) = 1.8371, vj = 5r s and a B = 0.05 as 
the fiducial value, hence, the free parameters are now M 
and £. 

The flux advection rate $ is determined from equa- 
tion (29) by specifying $ ou t. We assume that at w = 
ro ou t = 1000r s , £ out = Ik, T out = TWiai and /? out = 10, 
where T^iai = (/ic 2 /33?)(n7 out /r s ) _1 is the virial tem- 
perature at w = zu out . Since we fixed T, /?, and I at 
w — Tu ou t, equations (14) ,(15) and (18) indicate that 
W out cx £ ou t, H ont ~ constant and W out oc M, respec- 
tively. We find E ou t oc M, thus, equation (16) indi- 
cates that u roout ~ constant. Therefore, equation (28) 
gives B (n7 out ) oc Sou// 2 oc M 1 / 2 and equation (27) gives 
^outocM 1 / 2 . 

3.2. Numerical Results 

In figure 2, we show the thermal equilibrium curves at 
zu = 5r s for optically thin accretion disks. The equilibrium 
curves are plotted on the £ — M, £ — T, £ — P and £ — r c ff 
plane for £ = 0, 0.4 and 0.8. In the top panel of figure 2, 
M is normalized by the Eddington accretion rate MEdd = 
4TrGM/(r} e K es c) = 1.64 x lO 19 (O.l/r; c )(M/lOAf ) g s" 1 , 
where the energy conversion efficiency r) c and the elec- 
tron scattering opacity k os are taken to be r] e = 0.1 and 
k cs = 0.34 cm 2 g _1 , respectively. The equilibrium curves 
consist of three branches; (1) gas pressure supported, ad- 
vectively cooled branch, which is thermally stable, i.e. 
ADAF branch, (2) gas pressure supported, radiatively 
cooled branch, which is thermally unstable (SLE branch 
obtained by Shapiro, Lightman, Eardley 1976), and (3) 
magnetic pressure supported, radiatively cooled branch, 
which is thermally stable. We call the last branch "low-/? 
branch" . We find that the low-/3 branch exists even when 
M > O.lAfEdd, which corresponds to L > O.lLsdd- 

It may be worth noting that on the low-/? branch, M oc 
£ and T oc £~ 2 . This dependence can be explained as 
follows. When the magnetic pressure is dominant, W ~ 
B 2 H oc (v^B H) 2 /(v^ 2 H) oc &/{v^ 2 H) oc M/(v^ 2 H). 
We find w ro oc H^ 1 / 2 because equation (18) gives W oc M 
. Since equation (15) gives H oc {W/H) 1 / 2 oc (M/S) 1 / 2 



and equation (16) gives u ro oc M/£, we obtain M oc £ 
and H ~ constant. When the heating balances with the 
radiative cooling ( Q + = Q~), since Q + oc W and Q~ oc 
£ 2 T 1 / 2 / J ff, we find W oc M oc £ 2 T 1 / 2 /H. Using M oc £ 
and H ~ constant, we obtain T oc £~ 2 . 

The bottom panel of figure 2 shows the effective 
optical depth r c ff = ^K cs Kff £/2 where Kg = 6.4 x 
10 22 /?o(n7)T -3 - 5 cm 2 g -1 is the opacity of free-free absorp- 
tion. Note that, when £ < 0.3, the effective optical depth 
exceeds unity when M > O.lMEdd, thus, optically thin ap- 
proximation is no longer valid. We will extend the low-/? 
solution to optically thick regime in subsequent papers. 

One may wonder why T and (3 increase with £ when 
£ is fixed in the low-/? branch in figure 2. When £ 
is fixed, since W ~ {B /%ir)H oc B^W 1 ' 2 on the low-/? 
branch, Q + oc Vl^oc Bq 4 increases with £. Thus, the equilib- 
rium temperature increases according to W oc Q + = Q~ oc 
T x l 2 jH oc T 1 / 2 /W 1 / 2 as T oc W 3 oc Q+ 3 . The plasma (5 
on the low-/? branch is given by p ~ £T/W oc T 2 / 3 when 
£ is fixed. Thus, /? increases with £. 

Figure 3 shows the thermal equilibrium curves for «b = 
0.01. When a-Q is smaller, the transition from the ADAF 
branch to the low-/? branch takes place at smaller £ be- 
cause heating rate is smaller. The bottom panel of figure 
3 indicates that as the accretion rate increases, the disk 
becomes optically thick before the mass accretion rate at- 
tains M ~ O.lMEdd when £ < 0.5. 

In figure 4, we show the radial dependence of the M CY it 
above which the ADAF branch disappears for £ = 0.4 and 
a B = 0.05 (solid), 0.03 (dashed), and 0.01 (dotted). The 
critical accretion rate M cr i t increases inward. This in- 
dicates that as M increases the transition from ADAF 
branch to low-/? branch takes place first at large radius 
and propagates inward. 

4. Global Model 

We numerically solved equations (16)-(19) inward start- 
ing from the outer boundary at vo = vu ou t by specifying 
Lut, Tout, Pout, and four parameters a B , M, £, and £ in . 
By adjusting the parameter £; n , we obtained a global tran- 
sonic solution which smoothly passes the sonic point (e.g., 
Nakamura et al. 1997). The free parameters are now ae, 
M, and £. We adopt a B = 0.05. 

We locate the outer boundary at n7 ou t = 1000r s , and 
imposed the boundary conditions ^ ou t = 0.45^k, T out = 
^virial and p ou t = 10. We adopt backward-Euler method 
as an integration method (see Matsumoto et al. 1984). 

In figure 5, we show the radial structure of the disks 
when £ = 1 and M /M Ed d = 0.0002 (solid), 0.2446 (dashed), 
0.3181 (dashed-dotted), and 0.4094 (dotted). When M 
is small, since £ is small in the whole disk, the radia- 
tive cooling is inefficient. Thus the gas does not cool 
down. The heating is balanced mainly with the advective 
cooling. Furthermore, the flow is sub-Kcplcrian. When 
M is large, the disk can be divided into four regions; 
the outer boundary region w ^ 400r s , the outer transi- 
tion region 400r s ^ w ^ 150r s , the inner low-/? region 
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150r s > w > 50r s and the inner ADAF region w < 50r s 
when M/M-Edd = 0.3181 . Near the outer boundary, 
the accretion flow is advection dominated; Q a d v > Q~ ■ 
However, in the outer transition region, the radiative cool- 
ing overcomes the heating as E increases. Below this ra- 
dius T decreases inward, and /3 decreases, that is, the mag- 
netic pressure increases inward due to the magnetic flux 
conservation. At w ~ 100r s , T ceases to decrease inward 
because the heating balances with the radiative cooling. 
In the inner ADAF region, the advective cooling Q a dv in- 
creases as T increases, and balances with the heating Q + . 
The radiative cooling decreases in the innermost region as 
E decreases. 

To check the validity of the optically thin approxima- 
tion, we show the radial distribution of t c s in panel (f) of 
figure 5. In these models, T e s is less than unity. 

In the top panel of figure 6, we plotted the E — M re- 
lation at w = 5r s obtained by global transonic solutions. 
We also computed the luminosity L by integrating the 
radiative cooling all over the disk. In the bottom panel 
of figure 6, we show the dependence of L, which is nor- 
malized by the Eddington luminosity LecM = Vc^EddC 2 = 
1.47 x 10 39 (A//10M Q ) erg s _1 , on the surface density at 
vj = 5r s . The optically thin low-/3 solution extends to 
L ~ O.liEdd- Note that when E ~ 10 g cm~ 2 at w ~ 5r s 
the disk luminosity is only 1% of the Eddington luminosity 
even though the mass accretion rate attains M ^ O.lAfEdd- 
The disk is undcrluminous because the innermost region 
of the disk is still advection dominated, that is, radia- 
tively inefficient. The dotted curves in figure 6 (£ = 0.6) 
show that when 10 g cm~ 2 < E < 100 g cm~ 2 , L in- 
creases although M is nearly constant. It indicates that 
the disk luminosity increases as the area of the radiative 
cooling dominated region increases. The disk luminosity 
can exceed 10% of the Eddington luminosity only when 
M > M Edd - 

5. Summary and Discussion 

We obtained steady transonic solutions of optically 
thin black hole accretion flows by including magnetic 
fields. We theoretically confirmed the simulation result 
by Machida, Nakamura, Matsumoto (2006) that an opti- 
cally thin hot disk undergoes a transition to an optically 
thin magnetic pressure dominated disk when the accre- 
tion rate exceeds a threshold. In addition to the ADAF 
branch and the SLE branch known for optically thin disks, 
we found a low-/? branch in which the magnetic heating 
balances with the radiative cooling. 

In section 3, we obtained local solutions for single tem- 
perature disks dominated by the toroidal magnetic fields. 
As a consequence, we obtained three thermal equilibrium 
branches, ADAF, SLE and low-/? branches. 

Based on the results of global three-dimensional MHD 
simulations by Machida, Nakamura, Matsumoto (2006), 
we assumed that the Maxwell stress is proportional to 
the total pressure, (B vj B v /Att) = — aePtot- Here «b is a 
parameter corresponding to the a-parameter in conven- 
tional theory of accretion disks. Note that in magnet- 



ically dominated disks, since ptot ~ Pmag, the Maxwell 
stress is proportional to the magnetic pressure. In gas 
pressure dominated disk, since /3 ~ 10 inside the disk ex- 
cept in the plunging region of black hole accretion disks, 
Pmag ~ O.lptot- Thus our assumption is consistent with 
the results of global three-dimensional MHD simulations 
by Hawley, Krolik (2001, 2002), in which they showed that 
Maxwell stress is proportional to magnetic pressure (see 
also Pessah et al. 2006). 

In this paper, we assumed that «b does not depend on 
P and is uniform for simplicity. Let us discuss the depen- 
dence of «b on the strength of azimuthal magnetic fields. 
Pessah, Psaltis (2005) studied the stability of accretion 
disks with superthermal azimuthal fields and showed that 
MPJ is suppressed when the Alfven speed inside the disk 
va exceeds the geometrical mean of the sound speed c s 
and the Keplerian rotation speed «k (i.e. «a > \/c s vk)- 
Thus, the magnetic turbulence will be suppressed when 
this criterion is satisfied. In the three-dimensional MHD 
simulation reported by Machida, Nakamura, Matsumoto 
(2006), fluctuating magnetic field decreases as (3 decreases 
to /? ~ 0.1 (see Fig. 5b in Machida et al. 2006) because 
MRI is suppressed. Even in such disks, however, the an- 
gular momentum is transported by the Maxwell stress of 
mean magnetic fields, and as remained ae ~ 0.05 — 0.1 
(Machida, Nakamura, Matsumoto 2006). We need to 
carry out simulations for longer period to study the depen- 
dence of as on (3. We should note that «b increases in the 
innermost plunging region of the disk (e.g., Hawley, Krolik 
2002; Machida et al. 2004). Thus, this dependence may 
affect on the global solutions passing through the sonic 
point in the plunging region. However, its effect will be 
small for local solutions outside the plunging region. We 
should also note that we neglected the low-/3 corona above 
accretion disks. It will be our future work to take into ac- 
count the effects of such corona and outflows emerging 
from the disk. 

The heating rate inside the disk originates from the 
thermalization of magnetic energy via the magnetic recon- 
nection and is expressed as q + = (B^B v /47r)w(dn/dvD) = 
—Q!BPt tTZ>(dQ I 'dw) . This expression is based on the three- 
dimensional MHD simulations both in gas pressure dom- 
inant disks (Hirose et al. 2006) and in magnetic pres- 
sure dominant disks (Machida et al. 2006). If MRI is 
suppressed when va > y/c s VK as Pessah, Psaltis (2005) 
suggested, the dissipative heating rate Q + may decrease. 
Therefore, the disk temperature will become smaller than 
that shown in figure 3. 

As we showed in section 3.2, M oc E and T cx E~ 2 on 
the low-/3 branch. However, this dependence is differ- 
ent from that given by Machida, Nakamura, Matsumoto 
(2006) in which they derived M oc E 1 / 3 and T oc E~ 4 . 
This difference comes from the fact that in our model 
4> <x VmBH oc M 1 / 2 (that is, flux advection rate increases 
as mass accretion rate increases). Meanwhile Machida, 
Nakamura, Matsumoto (2006) assumed BH ~ constant. 

The critical accretion rate for the ADAF solution M cr ;t 
increases as decreasing the radius when w > 10r s , and has 
a maximum around w ~ 10r s . This fact indicates that in 
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the region vj > 10r s , the low-/3 region propagates inward 
as M increases. When M exceeds the maximum of Merit, 
whole disk will become low (3. 

According to our local solutions, optically thin low- 
(3 branches can exist above the critical accretion rate 
Af cr ; t above which the ADAF branch disappears (M cr ; t ~ 
O.lMEdd when «b = 0.05). The maximum luminosity of 
the disk staying in the ADAF branch is only 10 % of the 
Eddington luminosity. To explain the observed X-ray hard 
state with luminosity as high as 20% of the Eddington lu- 
minosity, the disk should be kept in optically thin state. 
Such disks can exist when £ > 0.5 in our models, where C 
is the parameter which parameterizes the radial variation 
of the advection rate of the toroidal magnetic flux. 

We also obtained steady global transonic solutions in- 
cluding magnetic fields and obtained E — M relation. 
When C ~ 0, we could not obtain solutions above certain 
critical accretion rate. The reason simply comes from the 
fact that we chose the ADAF-like outer boundary condi- 
tion. This means that there exists no solution connect- 
ing from such outer boundary to inner region with weak 
magnetic field. When £ ~ 1, we can obtain solutions for 
larger M, These solutions consist of four regions, the 
outer boundary region, the outer transition region, the 
inner low-/? region and the inner ADAF region. 

In this paper, we assumed single temperature plasma 
and included only the bremsstrahlung cooling. In mag- 
netized accretion disks, electrons will also be cooled by 
synchrotron radiation and inverse Compton effects. These 
processes will increase the luminosity of the disk and de- 
crease the critical accretion rate for the transition to the 
low-/? disk. In optically thin, hot accretion disks, since 
the gas temperature becomes so high while electrons are 
cooled by radiation, the electron temperature can become 
lower than the ion temperature. Such two-temperature 
accretion disk models were studied by Shapiro, Lightman, 
Eardlcy (1976) and Ichimaru (1977). Nakamura ct al. 
(1997) obtained transonic solutions of two temperature 
accretion disks including synchrotron cooling and inverse 
Compton effects. However, plasma (3 was assumed to be 
uniform in their paper. In subsequent papers, we would 
like to obtain solutions of magnetized two temperature 
disks subjecting to the cooling instability. 
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Fig. 1. A schematic picture of a magnetic field line inside 
the accretion disk, (a) Axisymmetric disk threaded by mean 
azimuthal magnetic fields, (b) Turbulent disk with mean az- 
imuthal magnetic fields and fluctuating magnetic fields. The 
azimuthal average of the radial component of the magnetic 
field, {B m ), is small because positive and negative B ra cancel 
out. On the other hand, the azimuthal average of the product 
of the radial and toroidal component, (B^Btp), has a large 
negative value because B-^B^ does not change sign when 
magnetic fields are deformed by nonlinear growth of MRI. 
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Fig. 2. Local thermal equilibrium curves of accretion disks 
at = 5r B on the £ — M, £ — T, £ — f3, S — r c ff plane for 
optically thin accretion disks with M = IOMq, qb = 0.05 and 
f = 0, 0.4, 0.8. SLE indicates the cooling dominated branch by 
Shapiro, Lightman, Eardley (1976). MEdd = ^GM/ (r] e K es c) 
is the Eddington accretion rate where the energy 
conversion efficiency r] c is taken to be r] c = 0.1. 
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Fig. 5. The radial structure of the disk when f = 1 
and M/M Edd = 0.0002 (solid), 0.2446 (dashed), 
0.3181(dashed-dotted), and 0.4094 (dotted) with 
«B = 0.05. Thick solid curve in panel (e) shows 
the Keplerian angular momentum distribution. 
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Fig. 6. (a) £ — M relation at 5r s obtained by global solu- 
tions with Ob = 0.05 when f = 1 (solid), 0.8 (dashed), and 0.6 
(dotted), (b) The relation between £ and luminosity L/Ledd- 



